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SUMMARY 

To explore the feasibility of using buckled columns to provide a soft 
support system for simulating a free-free boundary condition in dynamic 
testing, the nonlinear load-shortening behavior of initially imperfect, 
eccentrically loaded slender columns is analyzed. Load-shortening curves 
are obtained for various combinations of load eccentricity and uniform 
initial curvature and are compared, for reference purposes, with the 
limiting case of the classical elastica. Results for numerous combinations of 
initial curvature and load eccentricity show that, over a wide range of 
shortening, an axially loaded slender column exhibits load-deflection 
compliance which is of the same order as that of a straight but otherwise 
identical cantilever beam under lateral tip loading. 
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LIST OF SYMBOLS 

area of cross section 

constant of integration 

Young's modulus 

cross-sectional moment of inertia 

length of column 

moment 

applied vertical force 

radius of curvature 

moment arm for applied vertical force 

force-stiffness parameter 

parameter 

normalized vertical force 
parameter 

arc-length coordinate 
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t integration variable 

x vertical coordinate 

y horizontal coordinate 

7 column tip rotation 

A column shortening 

0 column cross-section rotation 

k column curvature 

Subscripts 0 and e denote initial state and Euler buckling value, 
respectively. 


INTRODUCTION 

In many structural dynamics tests it is desirable to minimize the 
interaction between a structure and its support system. One frequently 
used scheme involves suspending the structure on long, flexible cords 
when, for example, free vibration frequencies are of interest. Such an 
arrangement is most satisfactory when the suspension cords are long 
enough to minimize side loads and the natural frequencies of the 
suspension system are separated sufficiently from those of the test article. 

Another possible approach is to mount the structure on 
compression springs, if the springs can be made to have sufficiently 
low stiffness while supporting the weight of the test article. One device 
which might be able to satisfy the combined flexibility and load-bearing 
requirements, and simultaneously avoid the need for spacious overhead 
support capability, is the buckled column. After buckling, the column 
continues to support its buckling load, while over a considerable deflection 
range the slope of its load-shortening curve (i.e., its tangent stiffness) may 
be low enough to satisfy test requirements. 

To investigate the feasibility of such an approach, the shortening due to 
the large lateral deflections of an eccentrically loaded, slender elastic 
column with initial curvature is examined. Both the initial curvature and 
the load eccentricity serve to preclude buckling instability which would be 
troublesome in a support system. For reference, the load-shortening 
curves for various combinations of initial imperfection and load 
eccentricity are compared with results from reference 1 for the classical 
elastica problem. 



ANALYSIS 
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Governing Equations 

A sketch of a cantilevered, initially curved, inextensional column, before 
and during eccentric loading at the free end, is shown in figure 1. For 
simplicity, the initial imperfection shape is assumed to be a circular arc of 

radius R 0 . Hence, the initial tip rotation Yo is given by Yo=j^ ,where L is 

the length of the column. The vertical force P ia applied at a normal 
distance a from the column centroid, and so as to cause tip rotation of the 
same sense as the initial curvature. During deformation, the moment arm 
of length a is assumed to rotate with the end of the column. Equating the 
internal moment at a distance s measured along the column centerline 
from the tip, to the moment due to the applied force P yields 

-EI(K-K^)=P(y+acosy) /j 


where the initial curvature is Ko = -^-,Yis the angle of tip rotation relative 

Ko 

to the vertical, y is the lateral displacement relative to the tip, k is the total 
curvature, and El is the column bending stiffness. 


In terms of 0, the angle of rotation of the tangent at s, 
equation (1) can be written as 



Thus, 


dfi- = - J- - kfy + a cos y) 

QS Rq 


where k 2 = ^- . The boundary conditions are 
El 


and 


0(0) = y (and y = 0) 
0(L) = 0 


( 2 ) 

(3) 

(4) 


Differentiation of equation (2) gives 

d^6 = . k ^y 


ds 2 


ds 


or 



= - k^inQ 
ds 2 

Multiplication by dfi. and integration yields 

= 2k 2 cos 9 + C 

Applying boundary condition (3) and using equation (2) in 
equation (5) gives 


C = | + k^ cos yj 2 - 2k 2 cos y 

Thus, 

= 'V(r~ + k2a cos 2k2 ( cos Y - cos 0 ) 

where the negative sign is consistent with the fact that 0(s) is 
function of s. 


Employing the identity cos 0 = 1 -2sin 2 |, integrating over 
length of the column, and using boundary condition (4) yields 

r 


L = 


_d0_ 


Jo 


V (^- + k 2 a cosy) 2 + 4k 2 (sin 2 ^ 
V Kn 2 


sm% 


With the definition p = sin ^ , and the change of variable pt = sin 
(8) becomes 


4 

(5) 

( 6 ) 

(7) 

a decreasing 
the 

( 8 ) 

■, equation 



5 


where r = 


Yo + 


Tt^e q cos y\ 


. With the Euler buckling force for a 


4L I 


+ K^Tp^q 


perfect column denoted by p e =2t3E I, then k 2 =-^-q, where 


4L 2 


4L 


q = ■£- is the normalized applied force. Finally, equation (9) can be written 


as: 



dt 

VT ^ 2 fTr% 


( 10 ) 


which is the governing relation between the applied force and the initial 
and final tip rotations. 


To obtain the relation between applied force and the column 
shortening, which is defined herein as the vertical tip displacement in 
excess of that due to the initial curvature, an expression for vertical 
displacement is needed. (Note that in the present analysis the axial 

compressive compliance of the column has been neglected.) Since ^ = cos0, 

ds 

where x is measured downward from the tip of the column, then 


dx = cos0ds = cos0^d0 , so that use of equation (7) gives 
d0 


dx = - 




'sj + k^ cos yj 2 + 4k 2 jsin 2 ^ -sin 2 ®- 


01 ) 


With the definitions and change of variable employed earlier, integration 
of equation (11) over the length of the column yields 


X 


t 



V 1 - rk 2 
V 1 - pk 2 


dt -L 


( 12 ) 


where x t is the height of the deformed column, and equation (10) has 
been used to simplify the result. The column shortening due solely to 
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* * • * sin y 

initial curvature is L-L — ; hence, the net shortening due to the vertical 

To 

tip force is A = L - x t - [l - ij fo— ° | or 

\ Yo / 


A- lf sin Yo 4L.f VE S jt 
L Yo KVq J o 


(13) 


To obtain the relation between applied force and net column shortening 
for specified values of initial curvature and force eccentricity, equations 
(10) and (13) must be solved simultaneously. 


Computations 


A numerical approach was used in solving equations (10) and (13). 
Equation (10) was rewritten as 


f(Y.Yo>a>q) = 


2v( dL 

Vl-pWl-rV 


- 1 =0 


(14) 


which is the form required by routines for solving nonlinear equations. 
Before this equation could be solved, however, a numerical problem in 
evaluating the integrals in equations (13) and (14) required resolution. 
When r is very close to one, both integrals are nearly singular, giving rise 
to numerical integration inaccuracy. To circumvent this problem a change 
of variable was made, and asymptotic expansions were developed to 
evaluate the integrals on the subinterval 0.99 < t< 1.0. The integrals and 
their replacements are 

r i 


Ii = 


_dL 


Vi -pVfT- r 2t 2 


= ^(Jll + Jli 


where 


f 

Jn = 
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and 


where 


J 12= — 


VI 




+ : 


21 5 _LP1 9P 4 , 5P 6 
7V128 32 r 2 8f4 2 ' 


^V(l - T) 7 

r 6 / 


:=[ ^^dt = l{J 21 +Jd 
Jo Vl-rV 1 


Jai= 




vr^2 


du 


vri2 


jii= — 

22 VI 


! V 


and where, for the present computations, T = 0.99. 


For prescribed values of a, y 0 and y, equation (14) can be solved for the 

normalized force q. The integrals J n and J 2 i were evaluated using a 
standard mathematical library routine based on Simpson's rule. Equation 
(14) was solved by use of an available mathematical library routine which 
first bounds the root and then uses interpolation to obtain a converged 
value. An error criterion for f of 10‘ 4 was used in equation (14). After 

equation (14) was solved, A- was calculated from equation (13), again using 

the routine based on Simpson's rule. A complete load-shortening curve 
was generated by incrementing y and repeating the solution process for 
numerous prescribed values of a and y 0 . 


RESULTS AND DISCUSSION 

The results are presented in figure 2 as plots of q(4) for various 
combinations of initial curvature and load eccentricity. For reference, the 
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elastica solution, which the present solution approaches from below as load 
eccentricity and initial curvature vanish, is also shown in each figure. 

Figure 2(a) shows that for zero load eccentricity, the present results 
are barely distinguishable from the elastica solution when initial curvature 
is very small (Yo= 0.0001). As initial curvature is increased, the deviation 
from the elastica solution becomes more pronounced, particularly at small 

values of column shortening. However, for 4 greater than about 0.05, all of 

L 

the load-shortening curves are essentially parallel to the elastica curve, 
indicating that the tangent stiffness of all the "buckled" columns is about 
equal to that of the elastica. 

For progressively larger dimensionless load eccentricity, A , the relative 

importance of initial curvature is seen to diminish. In figure 2(f), for 
example, note the small range in the results for initial curvature values 
from 0.0001 to 0.03. Hence, a value of 0.05 for dimensionless load 
eccentricity would ensure smooth nonlinear load-shortening behavior, with 
little chance of troublesome column instability. However, the slope of the 

curve is somewhat greater than that of the elastica in the 4-range of 

interest, say, 0.1 to 0.3, indicating a tangent stiffness which is above the 
minimum achievable. 

Figures 2(b) through 2(e) are for values of load eccentricity in the 
range 0.002 < ^- < 0.02. As expected, all of these results lie in the envelope 

defined by the curves for load eccentricity, 4 , of zero and 0.05, and show 

in the 4-range of interest a gradual but discernable trend toward higher 

slope (thus, higher tangent stiffness) with increasing load eccentricity. 

Thus, the desire for the lowest practicable tangent stiffness would require 
that load eccentricity be kept reasonably small. 

The results also show that the effects of initial curvature and load 
eccentricity are essentially the same. Thus, the relatively soft 
"postbuckling" behavior of the nearly perfect slender column can be 
obtained by imposing reasonable values of either initial curvature or load 
eccentricity. In most situations it would appear that, of these two 



perturbing influences, load eccentricity would be easier to quantify and 
control. 

An appreciation of the magnitude of the tangent stiffness of the bent 
column can be gained by comparing the slope of the load-shortening curve 
for bending with that of the prebuckling load-shortening curve for the 
straight column under axial compression. For illustrative purposes, the 
two curves are combined in the sketch below to represent a continuous 
prebuckling and postbuckling load-shortening curve for a perfect column. 
The slope of the postbuckling curve is exaggerated for easier visualization. 



Note that tana = = £-[—] tanp « = tanp. 

\d A / 5 8 \L/ VdA/ c 


From the elastica solution (ref. 1), it can be shown (See Appendix) that 
the initial postbuckling slope of the load-shortening curve is, in 
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dimensionless terms, -^r = L which can be rewritten as (— ) = ^77. 

Ja_ 2 l dA/ b 8L 

where the subscript b denotes bending. (As can be seen in figure 2, this is 
a close approximation to the present results for a large range of column 
shortening.) The corresponding expression for prebuckling axial 

compression is f^-| = where the c denotes compression. The ratio of 

W c L 

the two tangent stiffness measures is 


(dP) 

_ ^2i _ 7i2/p\2 

WA l c 

where P is the radius of gyration of the column cross section. For slender 
columns, this ratio is very small. However, to put this result in 
perspective, note that the initial postbuckling slope is of the same order of 
magnitude as the slope of the load -deflection curve for the tip-loaded 
straight cantilever beam of length L, i.e., 

dP _ El 
dA 3L 3 

This means, of course, that a buckled-column support system would 
simulate free-free specimen support conditions about as well as would a 
system employing cantilever beams. Hence, although the slope of the load- 
shortening curve for the buckled column is very small in comparison with 
the perfect column's prebuckling slope, it is not negligible, particularly 
when compared with the initial postbuckling slope of the plot of load as a 
function of lateral tip deflection, which is zero,.,- 


CONCLUDING REMARKS 

In this paper, the load-shortening behavior of an eccentrically loaded, 
initially imperfect slender column has been analyzed. The results indicate 
that such columns can be designed to behave as relatively soft springs 



1 1 

over a wide range of column shortening, for reasonably broad ranges of 
initial curvature and load eccentricity. Both initial curvature and load 
eccentricity are seen to promote similar nonlinear behavior. Because load 
eccentricity is probably easier to impose with accuracy, it is likely to be 
the preferred parameter to specify in order to produce and control the soft 
"postbuckling" behavior of the slender column. If the tangent stiffness of a 
"buckled" column is low enough to satisfy test requirements, then a 
specimen support system that incorporates eccentrically loaded slender 
columns may present a useful alternative to a system requiring very long 
suspension cords. It should be noted, however, that the tangent stiffness 
of the deformed column is of the same order as that of an initially straight, 
tip-loaded cantilever beam and, thus, may not be low enough to satisfy 
some test requirements. 
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APPENDIX 


The Initial Postbuckling Slope of a Perfect Inextensional Column 

From reference 1, the two equations governing the load-shortening 
behavior of the elastica are, in dimensionless terms. 


f = ^-E(p)-l 
L kL 

(Al) 

kL = K(p) 

(A2) 


where K(p) and E(p) are complete elliptic integrals of the first and second 


kind, respectively, p=sin a , and k= ^-Vq". Letting \ = &-, (Al) and (A2) 

2 2 L 


become 


and 


hh=Me1 

K(p) 


1 2 

(A3) 


= K(p) 


Differentiating both equations with respect to 


K 


and 


dp = 

dX 




EdK. K dE 

dp dp, 


(A4) 


CA5) 


dq = iSdK dp 
d\ 71 dp dX 


Inserting (A5) into (A6) gives 


kAIK 

dq _ 2fq dp 

d\ n j$£lK - 

dp dp 

For load only slightly above the buckling load (q=l), p=sin 


SL 

2 


(A6) 


(A7) 

is very small. 


Hence, K(p) and E(p) have the valid power series 


and 


K( P ) = J 


1+ £- + -2-p 4 +. 

4 6 4 



(A8) 


(A 9) 


Substitution of (A8) and (A9) into (A7) yields 
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As q— »1+, p-»0; therefore, 
from above. 



1- as the buckling load is approached 
2 
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L 

Figure 2.- Load-shortening curves for an initially curved, eccentrically loaded column. 
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(b) = 0.002 

L 


Figure 2.- Continued. 
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Figure 2.- Continued. 




Figure 2.- Continued. 
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Figure 2.- Continued. 


(e) = 0.02 

L 




Figure 2.- Concluded. 


(f) A =0.05 
L 
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